Pheromone Static Routing Strategy for Complex Networks 
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Abstract 

In this paper, we adopt the concept of pheromone to generate a set of static paths that can 
reach the performance of global dynamic routing strategy [Phys. Rev. E 81, 016113(2010)]. In the 
test stage, pheromone is dropped to the nodes by packets forwarded by the global dynamic routing 
strategy. After that, static paths are generated according to the density of pheromone. The output 
paths can greatly improve traffic systems' overall capacity on different network structures, including 
scale-free networks, small-world networks and random graphs. Because the paths are static, the 
system needs much less computational resource than the global dynamic routing strategy. 
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I. INTRODUCTION 



Traffic on networked systems is important for many modern communication and trans- 
portation systems. Since the discovery of Small- World phenomenon and Scale- Free 
property [2] , it is widely proved that the topology and degree distribution of networks have 



profound effects on the processes taking place on these networks, including traffic flow [3Hl 1 ] . 

Routing is the control mechanism that directs each commuter in a networked system 
from its source to destination through a serial of switching nodes. A good routing strategy 



can improve network traffic efficiency without changing the underlying infrastructures 
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161 ]. Therefore, adopting efficient routing strategies are often preferred by the engineers. 
In modern communication and transportation systems, the shortest path protocol is often 
used. However, the shortest paths often lead to the collapse of hub nodes, which greatly 
reduces system performance. To avoid this problem, some routing strategies are proposed. 



One good example is the efficient routing strategy 



131 ] . The efficient path between any pair 



of nodes is defined as the path that the sum degree of nodes is minimum, denoted as: 



i 



Path = min'^^ k(x m )P , (1) 

m=l 

where k(x m ) is the degree of mth node along the path, I is the path length and (3 is a tunable 
parameter. When = 1, Eq.flT]) can reach a very high network capacity, which is almost ten 
times larger than shortest paths. Actually, the efficient routing strategy propels the packets 
to use the low-degree nodes, which composes the peripheral of the network. Therefore, the 
capacity of low-degree nodes can be utilized and so the system's capacity can be improved. 

The nodes' ability can not be fully reflected by the degree. The best routing strategy 
to fully utilize the nodes' ability might be accomplished by examine the queue length of 



nodes. Recently, Ling etal proposed a global dynamic routing strategy 



16J, in which the 



path between any pair of nodes is defined as the sum of nodes' queue length is minimum: 

l 

Path = rain [1 + n(x m )]; (2) 

m=l 

where n(x m ) is the queue length of node x m . This routing strategy selects the path along 
which the nodes have least traffic loads. Thus the network capacity can be further enhanced 
to almost double of the efficient routing strategy. As far as we know, the global dynamic 
routing strategy can achieve the highest performance for scale-free networks. Nevertheless, 



the shortcoming of this routing protocol lies in its high computational consumption since 
the paths are dynamical. This problem can be acute when the network size is large. An 
approach to reach the performance of global dynamic routing strategy with static paths is 
needed. 

In this brief report, we propose that by the help of "pheromone", this kind of static 
routing paths can be obtained. The concept of "pheromone" has been proved to be efficient 



to form a good local routing strategy 15|. Here it is demonstrated that "pheromone" can also 
help to generate routing paths that can reach high traffic capacities as the global dynamic 
routing strategy for scale-free, small-world and random network infrastructures. 

II. PHEROMONE STATIC ROUTING STRATEGY 



In the process of moving, ants wil 



release a special kind of secretion pheromone, which 



helps other ants to look for paths |l7l . Il8| . The pheromone released by the pioneer packets 
are important for the final paths to be established. Similarly, we suppose that some pioneer 
packets travel in the network by the global dynamic paths and drop an information of 
"pheromone" to the nodes in the first period. Therefore, the traffic flow characters of global 
dynamic routing strategy can be reserved. Then the static routing paths are generated by 
the "pheromone" information. Finally, the packets are forwarded with the output routing 
paths. The paths are static and the packets do not need to drop "pheromone" any more. 

The traffic model is described as follows. Initially, the pheromone density of each node 
is set to zero. The maximal queue length of each node is assumed to be unlimited and 
First-In-Fist-Out (FIFO) discipline is applied. The nodes' delivering ability is set to C = 1, 
i.e., each node can forward only one packet to its neighboring nodes in each step. 

(1) Test Stage: In this stage with preset length of T t step s, N pt test packets are forwarded 
in the system by the global dynamic routing strategy [161 ] . The routing paths are update 
in every step. Meanwhile, each test packet waiting in the queue will release one unit of 
pheromone to the current node for each step. That is, if a test packet waits for 5 steps 
in the queue before being delivered, the pheromone density of this node will increase by 
5 units. The unit pheromone value is set to 5p = 0.0001. When a test packet arrives at 
its destination, it remains in the system and is forwarded to another random destination. 
After T t steps, the nodes in the system will have different pheromone densities and the test 
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packets are all removed. 

(2) Path Generation Stage: The pheromone static path for any pair of nodes is generated 
that the sum pheromone density of nodes along the path is minimum: 

l 

Path = min p(x m ), (3) 

m=l 

where p(x m ) is the pheromone density of node x m . 

(3) Delivery Stage: In each step, R packets enter the system with randomly chosen sources 
and destinations. The packets are forwarded by the pheromone static paths. The paths do 
not need to be updated. Once a packet arrives at its destination, it is removed from the 



system. Following common practice, we adopt the order parameter J4]: 

v{R) = lip §^ (4) 

where AN p = N p (t + At) — N p (t), N p (t) is the number of packets within the network at 
time t, and (...) indicates the average over time windows of width At. With increasing 
packet generation rate R, there will be a critical value of R c which characterizes the phase 
transition from free flow (rj = 0) to congestion (77 > 0). In the free flow state, r) is around 
zero due to the balance of created and removed packets. When R exceeds R c , the packets 
accumulate in the network, and 77 becomes positive. The network's overall capacity can be 
measured by the maximal generating rate FL. R c can also be estimated analytically by the 
maximal node betweenness in the network [5j: 

R. - (5) 



■'max 



where betweenness (B) is defined as the number of paths passing through a node, and N is 
the network size. 

We note that this model is different from the ant colony algorithm 
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18} because it 



does not consider the evaporation effect of pheromone and the pheromone only appears in 
the first two stages. Actually, the concept of pheromone is borrowed only to record the 
traffic flow trace of global dynamic routing strategy. From the above model, one can also 
see that the model can be adjusted by changing the parameters of testing time length T t 
and test packet number N pt . In the following, we will show the results on three typical 



network structures: the Barabasi-Albert (BA) scale-free networks 



small-world networks 



19] and the Erdos-Renyi (ER) random graphs 



, the Newman- Watts 
2fJ. 
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FIG. 1: Pheromone density vs node degree k with test packet number N p t = 20000 and testing 
time Tt = 10000. Network parameters are size N = 1000 and average degree (k) = 4. 
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FIG. 2: (Color online.) Variation of system capacity R c with the length of test period Tt with 
Npt = 2000 (a) and Npt = 20000 (b). Each data point is obtained by carry out the model on one 
BA network realization with parameters N = 1000 and (k) = 4. 

III. SIMULATION AND ANALYTICAL RESULTS 

We firstly present the simulation and analytical results on BA scale- free networks j^J. In 
this model, starting from mo fully connected nodes, a new node with m links is added to 
the existing graph at each step according to the preferential attachment. The probability 
for the new node to be connected to an existing node i is proportional to the degree fcj of 
the node. 
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FIG. 3: (a) Maximal capacity R™ ax vs N pt with T t = 10000. (b) Saturation point T tc vs N pt for 
the system. Network parameters are TV = 1000 and (k) = 4. 
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FIG. 4: (a) Average queue length vs node degree k with the pheromone static routing strategy 
in free flow state, (b) The probability distribution of travelling time in free flow state (R = 60). 
Network parameters are TV" = 1000 and (k) = 4. 

Figure [1] shows a typical distribution of pheromone as function of node degree at the end 
of test stage. It approximately follows a piecewise power law with positive exponent. This 
indicates that the hub nodes are more burdened with the global dynamic routing strategy. 

Figure [2] shows the results of network traffic capacity measured by R c with changing T t . 
The analytical estimation of R c by Eq.Q are shown with Betweenness calculated by the 
pheromone static paths. For a given Npt, network capacity R c increases with T t until it 
reaches saturation at a critical point of T tc . The saturated network capacity R™ ax is affected 
by Npt- When N pt = 12000, the maximum network capacity is R™ ax = 68. This is very 
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FIG. 5: (Color online). Average travelling time for different R with network size N = 1000 and 
average degree (k) = 4. When R > R c , the average travelling time is calculated by averaging the 
travelling time of arriving packets during 1000 steps of delivery. 



close to that of the global dynamic routing strategy (R c = 71) [16J and is more than double 



of efficient routing strategy (R c = 33) [13J with the same network parameters. One can 
see that the analytical results are slightly larger than simulation results. This discrepancy 
might be due to the noise in the simulation. 

In Figj3](a), one can see that R™ ax increases with N pt and come to a saturate value of 68 
at Npt = 12000. In Figj3]^b), we also show the variation of critical point T tc for different Npt. 

Figure Eta) shows the distribution of average queue length as function of node degree k 
with the pheromone static paths. For both cases R = 60 < R c and R = 30 < R c , one can see 
that the traffic load distributes non-linearly with a maximum at k — 6. This indicates that 
the traffic load on hub nodes is alleviated. This behavior is different from the situations of 
shortest path and global dynamic routing strategy where the queue length follows a power 
law of n(k) ~ k 1 with exponent 7 > 0. Figure Hl^b) shows the probability distribution of 
packet's travelling time in free flow state (R = 60 < R c ) with the pheromone static paths. 
The travelling time is recorded as the time that the packet spend from source to destination. 



It approximately fo 
routing strategy 



16] . 



lows a Poisson distribution, which is similar to that of global dynamic 



Figure [5] shows the variation of packets' average travelling time over packet generating 
rate R for the efficient routing strategy and pheromone static paths, one can see that the 
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FIG. 6: Analytical estimation of network traffic capacity R c vs testing time Ti with N p t = 20000 
for: (a) Erdos-Renyi random networks of N = 1000 and (k) = 8; (b) Newman- Watts small-world 
networks of N = 1000 and (k) = 8. 

average travelling time is slightly larger than that with efficient routing strategy for R < 32. 
When R > 32, the system enters congestion for efficient routing strategy, thus the packet's 
travelling time increases rapidly. When R > 68, the system is also congested with pheromone 
static paths, thus the average travelling time also increases. 

Finally we present the results of network traffic capacity with Newman- Watts (NW) 



small-world networks 



19[ and Erdos-Renyi (ER) random graphs [20[. NW small- world 



networks are generated by randomly adding long-range links to a one- dimensional lattice 



with nearest and next-nearest interactions with periodic boundary conditions 19[. ER 



random graphs are generated by connecting couples of randomly selected nodes, prohibiting 



multiple connections [20J. We find the model can also achieve much higher system capacity 
on these networks. In FigJHl the analytical results of network traffic capacity measured by 
R c are shown for the two type of networks. The overall capacity can reach a maximum of 
R™ ax rs 234 for ER networks and R™ ax 222 for NW networks. For comparison, we find 
by simulation that R c ~ 50 with shortest path routing strategy and R c ~ 140 with efficient 
routing strategy for ER networks with the same parameters. For NW networks, R c ~ 85 
with shortest path routing strategy and R c ~ 155 with efficient routing strategy. 
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IV. SUMMARY 



In summary, a set of static routing paths is proposed for the traffic system on complex 
network infrastructures. Remarkably, the paths can almost reach the high performance of 
global dynamic routing strategy with a fixed routing table. It is demonstrated that the 
model can perform well on different network structures, include scale-free, small-world and 
random networks. For scale-free networks, the distribution of queue length and probability 
of travelling time are also investigated. This model can help the plan and design of routing 
protocols of modern communication and transportation networks. 
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